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ABSTRACT 

This  paper  develops  optimal  harvesting  policies  for  a  population  whose 
evolution  is  governed  by  a  single  autonomous  nonlinear  differential  equation. 
The  objective  functional  is  not  assumed  to  be  convex.  The  results  are  used  to 
discuss  optimal  policies  for  age-structured  populations  harvested  with  effort 
independent  of  age. 
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SIGNIFICANCE  AND  EXPLANATION 


This  paper  develops  optimal  policies  for  harvesting  biological 
populations.  The  age  structure  of  the  population  is  taken  into  account.  The 
policies  generally  involve  harvesting  at  maximal  effort,  not  harvesting,  and 
harvesting  to  keep  the  growth  rate  of  the  population  a  maximum. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 


summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 
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ON  OPTIMAL  HARVESTING  WITH  AN  APPLICATION  TO 
AGE-STRUCTURED  POPULATIONS 

*  ** 

Morton  E*  Gurtln  and  Lea  F.  Murphy 

1 .  Introduction. 

In  this  paper  we  develop  optimal-harvesting  policies  for  a  population  whose  evolution 
is  governed  by  a  single  autonomous  nonlinear  differential  equation.^  We  present  a 

2 

formulation  of  the  infinite-horizon  problem  which  does  not  Involve  a  discount  rate  and# 
using  straightforward  arguments,  are  able  to  establish  the  following  results: 

(a)  If  the  effort  is  unconstrained,  the  optimal  policy  is  to  reach  a  certain  value 
P*  of  total  population  as  quickly  as  possible,  and  then  to  hold  the  population  at  that 
value  for  all  subsequent  time.  As  one  would  expect,  P*  is  the  population  size  that 
maximizes  the  growth  rate  and  hence  corresponds  to  the  maximum  sustainable  yield. 

(b)  This  policy  is  also  optimal  when  the  effort  is  constrained  to  be  less  than  a 
constant  E,  as  long  as  E  is  sufficiently  large.  For  E  small,  however,  it  may  not  be 
possible  to  hold  the  population  P(t)  constant  at  P*.  In  this  case,  the  optimal  policy 
involves  a  "bang-bang*'  effort  function  and  the  total  population  approaches  a  limit  larger 
than  P* . 

We  prove  these  results  directly,  without  the  aid  of  classical  control  theory.  Indeed, 
the  classical  theory  is  not  well  suited  to  this  problem,  since  we  deal,  in  part,  with 
unbounded  controls,  and  since  our  formulation  of  the  infinite  horizon  problem  is 
nonstandard. 


Department,  of  Mathematics,  Carnegie-Mellon  University,  Pittsburgh,  PA  1S213. 

•  * 

Department  of  Mathematics,  Oregon  State  University,  Corvallis,  OR  97331. 

1 

This  problem  has  been  given  considerable  attention  in  the  literature  (cf.,  e.g.,  Clark  and 
Munro  f11,  Clark  (21,  Spence  and  Starrett  (31).  Our  results  (b)  for  and  E  small 

appear  to  be  new. 

2 

We  use  the  overtaking  criterion  of  optimality. 
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Science  Foundation  under  Grant  No.  MCS7fl-0193*i. 


We  use  our  results  to  discuss  optimal  policies  for  age- structured  populations 
harvested  with  effort  independent  of  age*  Using  a  simplified  form,  due  to  Coleman  and 


Simmes  [41,  of  a  general  nonlinear  model  introduced  by  Ourtin  and  MacCamy  (5),  we  are  able 
to  show  that  for  a  population  whose  age  distribution  is  initially  persistent^  -  an 
assumption  justified  for  a  population  which  has  been  evolving  over  a  long  period  of  time  - 
the  optimal-harvesting  problem  reduces  to  the  problem  discussed  above* 


2.  Age- independent  theory, 
a .  The  optimal-harvesting  problem* 

We  consider  a  species  whose  total  population 
by  a  differential  equation  of  the  form 


1.  The  cjrowth  rate  y. 


P(t),  when  not  harvested,  is  governed 

P  =  y(P) 

We  assume  that  the  growth  rate  y 
is  a  function  on 

with  a 

strict  global  maximum 

(1) 

at  >  0. 

We  assume  further  that,  on  (0,Pj,^), 
y  is  strlctly-positive^  and  has 
a  finite  number  of  critical 
points . 


1 

We  the  term  "p»=>rsistent  age  distribution”  for  what  in  the  literature  is  usually 

referred  to  as  a  “stable  age  distribution."  In  a  future  paper  [6]  we  will  investigate  the 
ronse:niencec  of  arbitrary  initial  conditions* 

This  asnumption  can  be  weakened  to  Include  "critical  depensatlon." 

'csiimo  that,  when  P,^  <  P„,  y  >  n  on 


Indeed,  we  need  only 


Suppose  that  at  t  =•  0  the  F>opulatlori  has  initial- value 

PfO)  -  Pg  >  0  ,  (2) 

and  that  subsequently  individuals  are  harvested  at  a  rate 

E(t)P(t) 

with  E  a  nonnegative  function  representing  the  (harvesting)  effort.  Then 

P  =  y(P)  -  EP  ,  (3) 

and  using  this  equation  and  (2)  the  yield 

/g  E(t)P(t)dt 

on  a  time  interval  !0,T)  can  be  written  as  a  functional  of  P: 

Y^{P)  -  Pq  -  P(t")  +  /g  y(P(t))dt  .  (4) 

Let  us  agree  to  use  the  term  path  for  a  ateictly-poaitive,  piecewise  continuous, 
right'continiious ,  ^  piecewise  function  P  on  [0,»).  Assume  that  we  are  given  a 

family  A  of  paths.  Then  an  optimal  path  relative  to  A  is  a  path  P*  6  A  with  the 
follow’';q  property:  given  any  PSA,  P  ^  P  , 

y^(p*)  >  y^(P)  (5) 

for  all  sufficiently  large  T  (that  is,  for  all  T  larger  than  some  T^  =  T^(P)).  Note 
that,  since  (5)  is  strict,  there  is  at  most  one  optimal  path. 

There  are  many  oossible  objective  (profit)  functions  with  which  one  can  work.  For 
converience,  we  have  chos*®  to  optimize  the  total  vleli.  Our  analysis  goes  through,  almost 
without  change,  when  a  cost  term  proportional  to 

1 

T.e.,  P(t)  =  P(t  ).  This  assumption  removes  ambiguity  at  points  of  discontinuity. 

p 

A  thorugh  discu.ssion  of  objective  functions  aoDropiiate  to  harvesting  problems  is  given  by 

Clark  f21. 
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/o  =<t)dt 


is  subtracted  from  <4). 

Remark .  The  inf inte-horlzon  problem  is  usnally  formulated  with  the  aid  of  a  discount 
”  (ft 

factor  e  (ic  >  0);  one  then  maximizes  the  total  yield 

/"  E(t)P(t)e''^dt  . 

This  formulation  renders  the  yield  in  the  immediate  future  more  important  than  that  of  the 
distant  future.  To  the  contrary,  our  procedure  models  situations  in  which  one  is  willing 
to  accept  a  sub-optimal  yield  initially,  in  order  to  eventually  produce  the  best  yield 
possible.  The  use  of  a  discount  factor  may  be  sound  from  an  economic  point  of  view,  but 
our  policy  is  certainly  more  appropriate  to  a  long  range  conservation  progr^un.  It  is 
interesting  to  note  that  discounting  generally  leads  to  a  smaller  size  for  the  ultimate 
population. 

b.  Optimal  harvesting  when  the  effort  is  allowed  to  be  unbounded. 

Our  first  step  will  be  to  define  the  class  A  of  paths.  To  begin  with,  we  require 
that  the  effort  be  nonnegatlvej  in  view  of  (3),  this  is  insured  by  the  constraint! 
for  any  t  >  0  at  which  P  is  differentiable, 

P(t)  <  y(P(t) )  . 

Next,  note  that  Y.j,(Pl  is  well  defined  for  all  paths  P,  even  those  that  suffer  jump 
discontinuities.  If  individuals  are  harvested,  rather  than  stocked,  the  population  after 
such  a  jump  must  be  less  than  the  population  before.  We  therefore  add  the  second 
constra int : 

(A^)  P(0)  ^  *by  t  >  0  at  which  P  jumps, 

P(t)  <  P(t")  . 

To  interpret  the  differential  equation  (3)  on  paths  with  jump  discontinuities,  we 


rewrite  it  in  the  form 


(  in  P> 


^  y(P)  -  E  . 


For  this  equation  to  make  sense  for  functions  with  jumps,  E  must  be  interpreted  as  a 
distribution.  Indeed,  if  tj^  (k  =  1,2, ...,n)  denote  the  times  of  discontinuity  of  P, 
then 


E(t) 


Eg(t)  +  5(t)in 


P(0) 


+  t  5(t-t^)Jln 
k  " 


P(t,  ) 


with  Eq  piecewise  continuous  and  6  the  delta  distribution.  Note  that,  with  E  defined 
in  this  manner,  (Aj)  follows  from  the  requirement  E  >  0. 

Since  y  has  a  strict  maximum  at  Pj,^,  we  would  expect  the  optimal  policy  to  involve 
reaching  P^^  as  quickly  as  possible  and  then  remaining  there  for  all  subsequent  time.  For 
Pg  >  P^  this  is  easily  accomplished  by  jumping  immediately  to  Pj^,  since  such  a  jump  is 
consistent  with  the  constraint  (Aj).  For  Pq  <  Pj^  a  jump  from  Pg  to  Pjj,  is  not 
possible)  as  we  shall  see,  here  the  optimal  way  to  reach  Pj^  is  to  refrain  from 
harvesting.  With  this  in  mind  we  introduce  the  following  definition:  the  f ree-growth 
curve  starting  at  Zg  e  (O.P^)  at  time  tg  >  0  is  the  solution  Z(t),  t  >  t^^,  of  the 
initial-value  problem 

Z  =  y(Z)  , 

Z(tp)  =  Z^  . 

Note  that,  since  y  >  0  on  Z  is  strictly-increasing  and  reaches  in  finite 

time* 

Theorem. ^  Relative  to 

A=  (paths  consistent  with  and  (^2^^ 

the  optimal  path  is  given  by 


1 

Cf.,  e.q.,  Clark  and  Munro  f1],  Clark  [2],  Chapter  2. 
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t  >  0 


(6) 


P*(t)  « 

for  ^  P  . »  and  by 

-  0  M  - ^ 


P  (t>  * 


Zit),  0  <  t  <  t 


P  ,  t  <  t  < 
.  M  M 


(7) 


for  P,-,  <  P^.  Here  Z  is  the  free«»qrowth  curve  starting  at  Pq  at  time  t  -  0,  while 
is  the  time  at  which  Z  reaches  P^* 

M  '■  I  —  —  ,  —  ■■■  ■■  —  M 


I'lqu'i:  -•  policies 


By  (3),  the  effort  E  required 
to  hold  P(t)  constant  at  Is 

given  by 

E  =  y(P„)/p„. 

Thus  for  the  path  (6), 

P 

E(t)  =  y(P„)/P„  +  6(t)in  —  , 
MM  P^ 

while  for  (7 ) , 


E{t) 


<  ':r  ’.  roof  of  the  theorem  begins  with 
,’\S5aeft  ton  1.  Let  P  0  A  be  optimal.  Then: 

I  i  )  tw  t  1  V  r  for  all  t  >0; 

(ii)  f'lr  all  sufficiently  large  t. 

Proof,  Assume  that  (i)  Is  not  true,  so  that  P  >  on  some  set  J  c  [0,“).  Define 

-  p  A  hv  )  -  min[P  ,Pft)).  Then,  since  G  <  P,  it  follows  from  (1)  that 

M 

Y  (T)  -  Y  (P)  =  '  "  y(P(t))]dt  +  P(t')  -  g(t')  >  0 

^  (0,T1  J  ” 


whenever  T  >  inf  J,  which  contradicts  our  hypothesis  that  P  be  optimal.  This  proves  (i). 
Before  beginning  the  proof  of  (ii)  note  that,  by  (1)  there  exists  an  L  e  (0,P„)  such 

M 

that 

P,  <  Pj,  Pj  e  (L,P„]  ->  y(P,)  <  yCPj)  .  (8) 

Our  first  step  in  establishing  (ii)  will  be  to  show  that 

lim  P(t)  -  P„  .  (9) 

M 


Assume,  to  the  contrary,  that  (^)  is  not  true*  Then,  by  (i), 

lim  inf  P(t)  <  , 

M 


and  there  exists  a  K  e  (L,Pj^)  and  an  increasing  sequence 


(t  }  with  t 
n  n 


and 


P(tj^)  <  K 

for  all  n.  Let  Z  denote  the  free-growth  curve  starting  at  K  at  t  =•  0,  let  k  denote 
the  time  at  which  Z  reaches  P^^,  and  let 


Z„(t)  -  Z(t-t„)  , 

so  that  Z.  is  the  free-growth  curve  starting  at  K  at  time  t_ . 


Further  assume,  without  loss  in  generality 
that  ^  ^n  * 

{t  )  does  not  have  this  property,  some 
n 

subsequence  does*)  Since 

P(t  )  <  Z  (t  )  , 

n  n  n 

Z  ==*  y(Z  ),  P  <  y(P)  , 
n  n 


Figure  3.  Boun(iinu  P  by  Z  . 

n 


and  since  the  jumps  in  P  are  downward, 
we  may  conclude  from  a  well  known  comparison 
theorem^  for  differential  inequalities  that 


1 

Cf.,  e.q.,  Hartman  f7],  Thm.  4.t,  p.  26. 
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(10) 


P  <  Z  on  (t  ,t  +k1 
n  n  n 

* 

Thus,  for  P  defined  by  (6)  or  (7)  and  sufficiently  large. 


[y(P*(t))  -  y(p(t))]dt  =  f  ty(P„)  -  y(P(t))ldt 
'  t  t  M 


>  L’'  fv(P„)  -  V(Z  (t))ldt  =  (y(P^)  -y(z(t))ldt  =  o  >  o 

t  n  *  I;  M 


where  we  have  used  (8):  hence 


F(T)  =  [yiP  (t))  -  y(P(t))ldt  +“  as  t  +  « 


Thus,  since  P  >  0, 

y^(P  )  -  y«(P)  =  F(T)  -►  P(T  )  -  P^  ♦  as  T  ®  , 

T  T  M 

which  again  contradicts  the  assumed  optimality  of  P*  Thus  P  satisfies 
Next,  by  (8)  and  (i),  there  exists  a  time  a  for  which 


Let  2  be  the  free-grnwth  curve 


starting  at  P(a^  at  tine  t  *  a, 


let  tc  denote  ♦•he  .•■ime  at  which  2 


reaches  P*,,  and  let  ^  A  be 
n 


i  defined  by 


Figuro  4-  'Construction  of  G. 


P( t )  ,  0  k  t  <  a 
G(  t)  =  (  Z(  t )  ,  a  C  h  <  < 


p^,  <  <  t  < 


-  Y^(P)  =  (y(n(t))  -  y(P(t))]dt  +  P(t")  -  P^ 


*'or  all  T  >  h-.  Fvirther,  the  argument  leading  to  (10)  here  tells  us  that 


P  <  Z  on  [aL,<] 


Suppose  that  (ii)  is  not  satisfied-  Then,  by  (R),  (11),  and  the  properties  of  G, 
y(G(t))  >  y(P(t))  for  all  t  €  [a,®)  , 

the  inequality  beinq  strict  on  an  open  set.  Thus,  since  P(t)  -*•  P  as  t  -»•  <»,  (  '  V  tiu 

M 

he  strictly  positive  for  T  sufficiently  larqe.  Since  this  cannot  be  so,  (ii)  must  be 
valid. 

In  view  of  Assertion  1,  it  suffices  to  establish  optimality  within  the  class 

A  a  {P  e  A  !  P  <  P,  on  P(t)  * 

MM  M 

for  all  sufficiently  large  tl 

*  * 

Assume  that  Pq  >  P^.  Let  P  be  given  by  (6)  and  choose  P  e  An,,  P  !<  P  .  Then 

P(T)  »  P*(T)  *  for  all  suf f icientl y  large  T,  and  for  all  such  T, 

V,,,(P*)  -  Y^(P)  =  /p  -  y(P(t)>ldt  >  0  ; 

* 

hence  P  is  optimal. 

To  complete  the  proof  we  must  establish  the  optimality  of  (7)  when 

Po  "  Pm  • 

With  this  in  mind,  we  introduce  the  following  notation.  Let  (A,C]  (h,P^),  lh*t  2 

denote  the  free-qrowth  curve  starting  at  A  at  t  -  0,  let  k'  denote  the  time  at  w'' i  c 
7  reaches  r,  and  let  [a,c^  (R,®)  with 


- 


c  >  a  +  1C 


(12) 


Further,  let  H  be  defined  on  la,c)  by 


H(t) 


( 


Z(t-a),  a  <  t  <  a  +  < 


a  +  <  <  t  <  c 


if  vtC)  >  y{h); 


H{t)  = 


a^t<c-ic 
Z(t-c+jcr),  c  -  <  t  <  c 


(13) 


!'i  .:nro  .  t}ptinuil  transitions.  Y(C)  <  y(A).  We  call  H  the 

optimal  transition^  from  A  to  C 
during  [a, cl j  <  is 

the  transition  time  of  H,  c  -  a  -  <  the  rest  time  of  H*  Finally,  a  path  P  0  is  a 

r  i  so  'inrina  fa.cf  if 

P{a)  <  P(t)  <  P(c)  for  all  t  e  (a,c)  .  (14) 

Assert  ion  2.  P  0  be  a  rise  during  (a, cl  Then: 

Iv)  There  is  an  .'ptimal  transition  H  between  P(a)  and  P(c>  during  (a, cl. 
(ii)  Tf  tl.e  rest  tine  of  H  vanishes,  then  P  *  H  on  (a,cl  . 

(in)  If  y  is  monotc'ne  on  lP(a)#P(c)l,  then 

y(H(t))>y(Pft))  for  all  t0fa,cl  , 
w'.th  i-.trirt  inefpiality  on  an  open  set  if  the  rest  time  of  H  is  nonzero, 
rtoii*  .  w*'  wi  \\  aivo  t  fu'  pr‘‘><*>f  onlv  for 

y ( p( c) )  >  v( P( a ) ) 

1  ♦  >  a  I  1  y  , 

P(t)  <  /.(t-a),  a<;t4a  +  <  (15) 


•  ’f,,  ypenc"  iind  5(tarr<'tt  f  1  ,  o.  303, 
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(cf.  (10)).  Thus  c  Is  consistent  with  (12)  and  the  optimal  transition  (13)  Is  well 
defined.  Assume  (for  the  remainder  of  this  paragraph)  that  the  transition  time  <  of  H 
equals  c  -  a.  Then  H(t)  «  z(t-a)  for  all  t  €  {a,c],  and  since 

H(c)  -  P(c)  , 

H  -  y(H),  P  <  y(P)  , 

It  follows'  that  P  >  H  on  ta.c);  this  inequality  and  (15)  yield  P  •=  H  on  (a,cl. 

To  prove  (ill)  let  y  be  monotone  and  hence  strictly  increasing  on  [P(a),P(c)]. 

(The  monotonicity  is  strict  since  y  has  at  most  a  finite  number  of  critical  points  in 
(0,P„).)  It  is  clear  from  (15)  and  (14)  that  P  <  H  on  la,c) ,  and  if  the  rest  time  of 
H  is  nonzero,  that  P{a+«)  <  Hta+ic)  =  P(c)e  This  yields  the  validity  of  (iii)s 

We  are  now  in  a  position  to  establish  the  optimality  within  and  hence  A  -  of 

the  path  P*  defined  by  (7).  Thus  choose  P  e  P  ^  P*,  and  let  be  the  least  time 

for  which 

P(t)  =  P„  for  all  t  >  .  (16) 

Since  all  jumps  in  P  are  downward,  P  must  take  on  all  values  in 

Cl''^2 . Cn  Po  =  ^1  <  =2  S  =  ‘'m 


I 


KX'tarf.’  6.  ixjcomposition  of  p. 


Hartman  (71,  Hemark  1,  p.  76, 


1 


denote  the  extreme  points  of 


n)  denote  the  last 


y  in  let  (k  =  1,2 

time  in  at  which  P  takes  on  the  value  Cj^  (note  that  let 

(k  =  2,3 . n)  denote  the  first  time  after  at  which  P  has  the  value  C^,  and 

define  =  0.  Then 

P  is  a  rise  during  t  i,  ,  T,  ,]  (k  =  1,2,...,n-1)  •  (17) 

k  k+1 

Consider  the  path  R  e  obtained  from  P  by  replacing  each  rise  (17)  by  the 

corresponding  optimal  transition.  Then,  since  y  is  monotone  on  each  [P(f^), 

we  may  conclude  from  Assertion  2  that 

y(R(t))  >  y(P(t))  for  all  t  >  0 

(Note  that  P  and  R  coincide  on  each  t  T,  ,  t,  1  • )  Thus 

k  k 

Y  (R)  >  Y  (P)  for  all  T  >  T„  .  (18) 

T  T  M 

We  shall  complete  the  proof  by  showing  that 

Y^(P*)  >  Y^(P)  for  all  T  >  T^  .  (18) 

*  • 

For  convenience,  we  consider  separately  the  two  cases:  R  =  P  ,  R  f  P  . 

Case  1  (R  =  P*).  Here 


(k  =•  1,2,...,n)  and  the  rest  time  - 
for  each  rise  except  the  last  -  vanishes  . 


(20) 


Thus,  by  (li)  of  Assertion  2,  R  and  P  coincide  outside  of  the  last  rise-interval 

[i  ,T  ).  The  rest  time  associated  with  this  last  rise  cannot  vanish,  because  if  it  did, 
n- 1  n 

* 

then  R  and  P  would  coincide  everywhere,  an  impossibility,  since  R  =  P  /  P. 

Consequently,  we  may  conclude  from  (iii)  of  Assertion  2  and  (8)  that 

Y™(P>  <  Y„(R)  (=  Y„(P  ))  for  all  T  >  T  . 

T  T  T  M 

Case  2  (R  ^  P*).  Let  T  >  T„  =  f  .  By  (7),  the  integral 
-  M  n 

y(P  (t) 
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* 

is  equal  to  ( T-Tj^ ) y ( )  plu'  the  inteqral  of  v(P  )  over  the  free-qrowth  segment  of 

P*.  Since  the  free-growth  equation  is  autonomous,  this  latter  integral  equals  the  integral 

of  y(R)  over  the  free-qrowth  portions  of  R.  Thus 

Y,j,(p'')  -  Y,y(R)  =  (T-T„)y(P„)  -  I,j,(R)  ,  (21) 

where  I,j,(R)  is  the  inteqral  of  y(R)  over  the  set  Q  consisting  of  the  rest  portions  of 

the  optimal  transitions  that  comprise  R  and  the  union  of  the  intervals 

( £  ),...,(t  ,T).  On  f),  y(R)  <  and  (as  we  shall  show  at  the  end  of  the  proof) 

2  2  n  M 

(1  contains  an  open  set  0  with 

y{R(t))  <  y(P„)  for  all  t  e  0  .  (22) 

Thus,  since  the  measure  of  Q.  is  T  -  Tj^, 

I,j,(R)  <  (T-T„)y(P„)  , 

and  (21)  yields  Y,j,(P*)  >  Y,j,(R).  In  view  of  (18),  this  inequality  implies  (19). 

To  verify  (22)  note  first  that  (20)  cannot  be  satisfied,  for  if  it  were,  then  R 

would  equal  P*.  Thus  at  least  one  interval  or  at  least  one  rest  interval  with 

R  P.,  has  nonzero  duration.  This  interval  clearly  contains  an  open  set  0  on  which  (22) 

is  satisfied.  (If  the  interval  In  question  is  (r  1#  the  existence  of  0  is  insured 

n  n 

by  the  equality  of  R  and  p  on  fact  that  is  the  smallest 

time  consistent  with  (16)). 

Remark  1.  The  construction  P  -  obtained  from  P  by  replacing  each  rvse  (between 
extrema  of  y)  by  its  corresponding  optimal  transition  -  can  he  defined  (in  an  obvious 
manner)  for  anv  function  P  with  domain  an  interval  fa,hl,  as  lono  as  P(a)  <  Pfb)  and 
P(t)  <  P(b)  for  all  t  ^  (a,b).  We  will  refer  to  this  construction  as  the  guas  i  ni  r.er 

nf  p  f  a  ,  b’  , 

c •  Bounded  effort. 

The  problem  becomes  more  coir.r.lex  when  the  effort  is  restricted  by  an  upper  bouni)  F. 

We  require  that  E(t)  v,  E  for  all  t,  or  equivalently, 

fA^)  P(t)  >  Y(P(t) )  -  EPft) 

A  path  whi^h  obeys  (A.j)  and  (Aj)  necessarilv  continuous.  We  therefore  Msmiss  ml 

re.strict  our  attention  to  paths  consistent  with  (A^),  (A-j),  and  the  initial  condition  (2). 


In  addition  to  assuming  that  y  Is  a  function  on  [0,«o)  with  a  finite  number  of 

critical  points,  we  assume 

lim  sup  y(P)  <  0  .  (23) 

P  ♦  09 


Then  there  is  a  largest  value  which 

y(Ps)  =  iPg  . 

Further,  by  (23),  y  restricted  to  ®  maximum;  we  assume  that  this  maximum 

occurs  at  exactly  one  point  P,.  e 

M  S 

yfP^  )  >  y(p)  for  all  p  >  p  ,  P  P„  •  (24) 

M  S  M 

Finally,  we  assume'  that  y  >  0  on  (0,Pg>.  We  use  the  notation 

S  =  yCPg),  M  =  y{P„)  ;  (25) 

then  M  >  S. 

In  addition  to  free  growth  curves,  the  optimal  path  involves  maximal-effort  curves; 

that  is,  solutions  X(t),  0  <  t  <  ■»,  of  the  equation 

X  =  y(X)  -  EX  .  (26) 

2 

We  again  appeal  to  a  comparison  theorem  to  note  that  if  P  satisfies  (A3)  and  Pttp)  • 
X(tQ),  then 

P  <  X  on  (0,tg),  P  >  X  on  .  (27) 

For  Pj,  <  Pj^,  we  will  use  the  following  notation: 

=  inf{p  >  Pp  !  y(p)  >  M)  , 

and  for  each  value  of  n(n  =  1,2,...)  for  which  the  underlying  set  is  nonempty, 

V  =  inf(p  >  U  (  y(p)  <  m)  , 
n  n 

U  =  inflp  >  V  I  y(p)  >  M} 
nF  I  n 


Then 


1 

: ;  i  :  i^r  *  i''n  '^.ly  be  wpakcne-i  as  Indicated  in  Footnote  2,  p.  2* 
Har*:man  [7],  Thm*  4.1  and  Remark  1,  p*  26. 
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(28) 


y  >  M  on  (U  ,V  )  , 

n  n 

and,  by  (24),  V  <  P  ,  for  all  n.  Since  y  )ia8  a  finite  number  of  extrema,  it  is  clear 
n  S 

that  there  is  a  finite  number  N  (say)  of  0_*s  and  the  same  number  of  V_'s. 

^  n  n 


We  consider  two  cases. 


Case  A,  P|^  is  not  a  local 


maximum  of  y.  Then,  by  (24), 


P„  =  Pg,  y  <  S  on  <Pg»“)  r 


^0  “l  ''l 


Figure  7.  The  function  y  for  Case  A. 


Theorem  A.  Relative  to 


p  and  Vy. 


A  »  (paths  consistent  with  (A^),  (A^),  and  (2)) 


the  optimal  path  is  P  ,  as  described  below. 


(a)  If  Pf,  <  P„,  then 


I'  0  <  t  <  u, 

I  X^(t),  u^  <  t  < 


Vi  ‘  <  "n 


X^(t),  <  t  < 


the  free-qrowth  curve  through  Pq  at  t  =  0, 

the  free-qrowth  curve  through  ^ 

the  time  at  which  2^.  (u^)  = 
n  -  n  n  n 


the  time  at  which  X^(v^)  - 
- - - n  n  n 
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(b)  ^  Pg  >  P  ,  then  P  la  the  maximal-effort  curve  through  Pg  ^  t  »  0. 


We  justify  the  by  noting  that, 

for  n  <  N,  y(p)  -  Ep  >  0  on  , 

so  that  X_  reaches  V„  in  a  finite 
n  n 

time.  The  optimal  paths  are  graphed  in 
Figure  9.  Before  proving  the  theorem  for 
Case  ft,  we  state  the  results  for 


U  VP 
N  NS 


Figure  8.  The  function  y  for  Case  B. 

Case  8.  Pj^  is  a  local  maximum  for  y  (Figure  8). 

Theorem  B.  Relative  to 

A*  (paths  consistent  with  (ft^),  (ftj),  and  (2)} 
the  optimal  path  is  P* ,  as  described  below 


(»  11  Pq  <  Pm' 


Z,(t),  0  <  t  <  U, 


X,(t),  u,  <  t  <  V, 


''n  "  "  ^  "n 


‘  ^ 


Here  Z^,  X^ ,  u^^  and  are  Interpreted  as  in  Theorem  ft,  is  the  free-growth  curve 

through  at  t  =  v^,  and  t^,  is  the  time  at  which  Z^^,(t^)  =■  Pn,. 


^0  »  Pm  ’  Pg’ 


X(t) ,  0  <  t  <  t 


P«'  ^ 


where  X  is  the  maxlmal~ef fort  curve  through  ^  t  *  0,  and  is  the  time  at 

which  X(t„)  *  P^,  (Since  X(«)  *  P_  <  t  <  ®. ) 

-  M  M  —  —  S  MM 

(c)  T_f  and  Pg  =  Pj^j,  then  P*  Is  the  maximal-effort  curve  through  P^ 


at  t  **  D. 
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Figure  9.  Optimal  paths. 

We  will  prove  the  theorem  only  for  Case  A(a)*  Since  the  technical  aspects  of  the 
proof  are  rather  complicated,  we  offer,  as  a  preliminary,  the  following  intuitive  argument 
in  support  of  the  optimality  of  the  path  (30). 

We  first  ask,  is  there  an  optimal  equilibrium  value  for  the  population?  By  (A^),  any 
equilibrium  value  p  must  satisfy  v(p)  <  Ep.  It  is  clear  from  Figure  8  that  the  "best” 
such  value  is  p  =  Pg.  We  therefore  expect,  and  will  indeed  show,  that  the  optimal  path 
asymptotically  approaches  Pg. 

We  next  ask,  what  Is  the  best  way  to  make  the  population  crow  from  Its  current  size  to 
the  optimal  size  P^?  We  know  that  the  optimal  asymptotic  yield  rate  is  S.  Thus  if  the 
current  yield  rate  is  less  than  S,  we  refrain  from  harvesting  lo  hasten  the  growth  of  the 
population  to  a  size  which  will  produce  a  higher  yield.  On  the  other  hand,  when  the 
current  yield  rate  is  higher  than  we  are  in  an  unusually  profitable  situation. 

Although  the  population  will  inevitably  grow  beyond  this  nrofitablo  stage,  we  can,  hv 

* 

harvesting  with  maximal  effort,  prolong  the  period  of  high  profit.  The  path  P  as 
defined  in  (30)  behaves  in  the  manner  just  described. 
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which  will  be  the  case  if 


Remark  2.  If  y  has  a  strict  global  maximum  at  >  Pg, 

E  is  sufficiently  large,  then  it  is  possible  to  harvest  so  that  P  is  held  constant  at 
Tor  this  case  Theorem  B  is  completely  analogous  to  our  theorem  for  unbounded  effort. 
For  Pjj  <  Pj^  the  optimal  path  is  a  free-growth  curve  followed  by  P(t)  5  P^;  for  Pg  >  Pj^ 
the  optimal  path  is  a  maximal-effort  curve  (which  replaces  the  jump  to  followed  by 

P(t)  ;  P^. 

We  now  begin  our  proof  of  Theorem  A(a>e  Thus  let 

Pq  <  • 

Assertion  3.  Let  PSA  be  optimal.  Then: 

( i)  lim  P(t)  *  Pgj 
t  -^0“ 

(ii)  P  is  nondecreasing. 

Proof .  Note  that  P(t)  <  P^,  0  <  t  <  ",  since  otherwise,  the  path  Q  e  A  defined  by 

Q(t)  =  min{P(t),Pg} 

would  be  better.  (For  convenience,  we  use  the  terminology:  Q  is  better  than  P  if 
Y^(C)  >  Y^(P)  for  all  sufficiently  large  T.)  Assume  first  that  ®  come 

time  tg.  Let  X  be  the  maximal-effort  curve  through  Pttg)  at  tg.  By  (26),  X(t)  •*  Pg 
as  t  -►  «:  thus,  since  X  <  P  <  Pg  on  ttg,"),  ( i)  is  satisfied.  Assume  next  that 

P(  t)  <  Uj,  for  all  t.  Lot 

Sg  =  {p  e  (0,0^1  I  y(p)  >  s  -  e)  .  (3i) 

Then  there  exists  an  e  >  0  such  that 

y(p)  >  Ep  for  all  p  6  .  (32) 

Moreover,  since  y  has  at  most  a  finite  number  of  critical  points,  Sg  is  the  union  of  a 
finite  number  of  compact  intervals.  By  (26),  (27),  and  (32),  if  P  enters  any  such 
interval,  it  remains  only  a  finite  time,  and  it  never  returns.  Thus  P(t)  lies  outside 
Sq  for  all  sufficiently  large  t,  and  (25)^  and  (31)  yield 

y(Pg)  -  y(P(t))  >  e  (33) 

for  all  such  t.  T,et 


-18- 


Q(t)  - 


0  «  t  < 

s 

t„  <  t  <  <• 

s 


where  z  is  the  free-qrowth  curve  etartlnq  at  Pp  at  t  ■  0,  while  tg  is  the  time  at 
which  z  reaches  Pg.  Then  a  trivial  computation,  based  on  (21),  shows  that 

lim  tY^(Q)  -  Y^(P)1  -  +<•  , 


so  that  Q  is  better  than  P.  This  contradicts  the  optimality  ot  t>  and  the  proof  of  (1) 
is  complete. 

♦  + 

(ii)  We  will  show  that  P(t  )  >  0  for  all  t*  Assume^  to  the  contrary,  that 

•  4-  — 

P(t^)  <0  at  some  time  t,|,  and  let  *  P{t^>.  Then  by  (A^),  y(P^)  <  EP^ ,  and, 
since  Y^P5)  =  ^  ^S'  addition, 

q  »  aup{p  e  [0,P^)  I  y(p)  =  Ep)  , 

T  =  inf{p  e  (P^.Pgl  I  Y<P)  ”  Ep} 

exist  (recall  y(0)  >  C),  and 


y(p)  <  Ep  for  all  p  6  (q,T)  , 

0  <  q  <  P^  <  T  <  Pg  . 


(34) 


(Cf.  Figure  10,  The  inequality  T  <  is  a  consequence  of  (28)  and  (34)^). 

Let  X  be  the  maximal-effort  curve  starting  at  P^  at  t  *  t^.  Since  the  maximal- 

effort  curve  through  q  remains  constant  at  q,  X  >  q  on  hence  (27)  yields  the 

conclusion 

P>qon(t^,a>),  (35) 

Let  t^  +  i  be  the  first  time  after  at  which  P(t^+i)  =  P^.  (Since  <  Pg, 

the  existence  of  2,  >  0  is  assured  by  (1).)  Then 

P  <  on  (t^,t^  +  )l)  ,  (36) 

By  (34)2,  (35),  and  (36),  q  <  P(t)  <  t  on  (t^,t^+i),  so  that,  by  (34)^, 

Y(P(t))  <  EP(t)  <  Et  =  y(T)  f-^r  all  t  e  (t^,t^  +  £)  . 

Thus  the  path  0  8  A,  pictured  in  Figure  11  and  defined  by 
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0(t) 


'  Q(t),  0  <  t  < 

P(t+t),  <  t  <  -  £ 

]  T,  -  £  <  t  < 


V  p(t) ,  <  t  <  ®  , 


where  t2  (>  t^+£)  is  the  first  time  after  t^  at  which  ^(^2)  *  T/  is  better  than 

•  + 

Hence  P(t^)  <0  is  not  possible.  This  completes  the  proof  of  Assertion  3. 


Figure  10  Figure  11 


The  foUowinq  definition  will  be  useful*  Let  f  be  a  continuous,  nondecreasino 
function  on  and  let  F  belong  to  the  range  of  f.  Then  f  stretched  at 

F  by  amount  i  >  0  is  the  function 
I  J  defined  by 


f  /  r  f<t),  0  ^  t  c  t^ 

F  -  ^ - ^  f(t)  =  \  ’'O 

I  f(t-£),  t„  +  £  <;  t  <  » 

y  0 

- 1 - 1 

0  0 

Figure  12.  f  stretched  at  F. 
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where  tQ  is  the  first  tine  at  which  *  F*  A  typical  stretched  function  is  shown  in 

Figure  12.  Of  course «  by  successive  applications  of  this  definition  we  may  stretch  a 
function  at  a  finite  number  of  points  in  its  range. 

* 

We  are  now  in  a  position  to  establish  the  optimality  of  the  path  P  defined  by 

* 

(30).  In  view  of  Assertion  3,  it  suffices  to  show  that  P  is  better  than  any  path 

P  e  A/  P  ^  P  /  consistent  with  (i)  and  (ii)  of  Assertion  3.  Let  P  be  such  a  path/  let 

u  (n  »  1,2,.. ./N)  denote  the  time  at  which  P  crosses  U_,  and  let  v  (n  =  1,2,...,N-1) 
n  n  n 

denote  the  time  at  which  P  crosses  (These  times  ate  unique?  Indeed,  if  t  denotes 

any  one  of  these  times,  then  P(t)  >  y(P(t))  *•  EP(t)  >  0.)  Let  P  be  P  stretched  at 

.  with  X  -  the  amont  of  the  stretch  at  V_.  Here  X  is  the  time 

it  takes  the  maximal-effort  curve  starting  at  to  reach  V^.  Finally,  let  u^ 

(n  »  1,2,«..,N)  be  the  time  at  which  P  crosses  U  ,  so  that  v  »  u  +  X 

n  n  n 

(n  »  1,2,...,N-1)  is  the  last  time  P  takes  on  the  value  V^.  Then,  since  P  ^  P  , 

P  7*  P*  I  (37) 

in  addition,  a  simple  analysis  shows  that 

lim  fy^(P)  -  Y^(P)1  =  0  . 

T-*«o 

It  therefore  suffices  to  show  that 

*  ♦ 

P  is  better  than  P  ,  (3H) 

♦ 

Consider  the  path  R  obtained  from  P  as  follows; 

(I)  On  [u  ,v  ]  replace  P  by  the  maximal-effort  curve  starting  at  at  t  =  u  . 

n  n  n 

(Note  that  P  is  defined  so  that  this  maximal-effort  curve  reaches  at  v  .) 

n  n 

(IT)  On  (v  .,u  1  replace  P  by  its  quasi-optimizer  (cf.  Remark  1).  Here  (I)  and 
n- 1  n 

-<*  -► 

(II)  hold  for  n  =  1,2,...,N  with  v^  =  0  rin  ?  v  =  «>, )  The  paths  P  and  R  are  shown 

”  N 

in  Figure  13.  Note  that  these  oaths  gener-^lly  do  not  belong  to  A. 


Figure  13.  The  curves  P  and  R. 


Assertion  4. 

lim  inf  (Y^(R)  -  Y^(P)J  >  0 
T  oo 


with  strict  tnpcuality  if  R  P. 


Proof «  Let  T  >  and  define 


T  =  U  fv  .u  J  =  D  fu  ,v  ] ,  j  ^  j  0  [0,T] 

,  n-1  n  ,  n  n  T 

1  r\~^ 


HA)  =  [y(R(t))  -  y(P{t))]dt  , 

A 


(3P 


so  that  (0,T]  =  I  u  and 


Y^(R)  -  Y^fP)  =  P(T)  -  R(T)  +  I(I)  + 


rfn  )  =  R(vi  ),  f?7)^  implies  P(T)  >  R(T)%  Thus  to  prove  Assertion  4  it  suffices 


to  show  that  (40)  and  (41)  below  are  valid: 

Id)  >0  with  strict  inequality  if  R  P  on  I  ;  (40) 

lim  inf  ?{J^)  >  0  with  strict  Inequality  if  R  P  on  J  .  (41) 

T  -*■  ® 

Conclusion  (40)  is  an  immediate  consequence  of  (ii)  and  (iii)  of  Assertion  2  and  the 
fact  that  R  P  on  I  only  if  at  least  one  of  the  optimal  transitions  associated  with 
R  has  a  nonzero  rest  time. 

We  now  prove  (41).  On  an  interval  tu  ,v  ) ,  R  is  a  maximal  effort  curve.  Because 

n  n 

♦  -►  — +  -V 

P  >  y(P)  -  EP  >  0  whenever  U  <  P  <  V  ,  P  must  be  strictly  increasing  on  [u  ,v  ] 

n  n  n  n 

until  the  time  at  which  P  reaches  Of  course,  if  n  ■  N,  P  may  never  reach  Vj^,  in 

which  case  P  is  strictly  increasing  on  (u^,“).  Let  C^,  <  Cj  <...<  »  V^, 

denote  the  extrema  of  y  on  so  that  K  is  odd,  and 

for  k  even,  y  ts  strictly  increasing  on  ' 


(42) 


strictly  decreasing  on 


Let  t.  be  the  time  at  which  R(t.  )  *  C.  .  Further,  let  Q  be  the  function  obtained  by 


stretching  P  at  each  for  odd 

k  <  K,  and  then  restricting  the  result 

to  [u  ,v  ] .  The  amounts  of  the  stretches 
n  n 

are  uniquely  determined  by  the  requirement 
Q(tj^)  =  R(t^)(=  C^)  for  k  even.  It 
is  clear  from  this  construction  and  (27) 
that,  for  k  even, 


r 

'k-l 


<  Q  ^  R  < 


on 


[t, 


<  R  <  Q  4  C 


k+1 


lnpauali^i«s  and  ^42)  im’" y  that 

y(R(t))  >  v(Q(t))  on  (u  ,v  1 

n  n 


(43) 


The  integral  of  Y(P)  over  the  times  at  which  P  is  strictly  increasing  is  equal  to  the 


analogous  Integral  for  y(Q).  Thus,  since  y(C.  )  >  y(C  ), 

K 

V 

!j'  ty(Q(t))  -  y{P(t))]dt  >  0  , 

u 

n 

and  we  conclude  from  (43)  that 

V 

fj'  (y(R(t))  -  y(P(t))]dt  >  0  , 

u 

n 

for  n  <  N.  Thus  to  show  that  lim  inf  I(J^)  >0  we  have  only  to  show  that 

lim  inf  [y<Q(t))  -  y(I^(t))ldt  >  0 

T  ♦  «  U 

n 

To  verify  (45)  let  ()<  odd)  be  the  length  of  the  stretch  in  Q  at  and  let  L 

denote  the  sum  of  the  Then  the  integral  in  (45)  is  equal  to 

T-L  ^  • 

odd  k 

The  integral  tends  to  zero  as  T  +  oo,  since  P(t)  ♦  Pg  as  t  ■*•  <»,  while  the  sum  is 
>  0,  since  Thus  (45)  holds. 

The  remainder  of  the  proof  of  (41),  that  lim  inf  I(J,j,)  >  0  when  B  P  on  J, 
safely  be  omitted.  This  completes  our  proof  of  Assertion  4. 

We  now  return  to  the  proof  of  (38)  and  consider  separately  two  cases. 

*  > 

Case  1  (R  *  P  ).  Here,  by  (37),  R  P,  and  (38)  follows  from  Assertion  4  with 
R  =  P  . 

♦  -w 

Case  2  (R  ^  P  ) .  Choose  T  >  u^  (>  u^)  and  consider  the  integral 
Ig  =  Jg  {y(P*(t))  -  v(R(t)))dt  . 


(44) 


(45) 


can 
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The  inteqral  of  y(R)  over  the  free-qrowth  portions  of  R  is  equal  to  the  analoqous 


* 

inteqral  for  P  .  Similarly,  the  inteqral  of  y(R)  over  the  maximal-effort  secjments  of 
R  on  [n,u^]  equals  its  counterpart  for  P  on  [O^u^],  Thus  ^0  ”  ^2  ”  ”  ^3'  where 

^  "  /I  ylRIt))'’*:,  =  /I  y(p*(t))dt 

“n  '^n 

and  Is  the  inteqral  of  y(R)  over  those  subintervals  of  I  (cf.  (39)^)  on  which  P 

is  constant*  The  measure  of  these  sub-intervals  Is 

^  \  '  “n  • 

Thus,  since  y(R(t))  <;  S  on  I,  we  have  the  bound 

<  tS  -  A 

with 


( 1 ) 

A  =  0 

when 

y(R(t) ) 

=  s 

on 

all  of  the  above  subintervals 

(2) 

A  >  0 

when 

Y(R(t) 1 

<  s 

on 

at  least  one  such  subinterval 

Finally, 

since 

R(t+T) 

=  P*{t) 

for 

t 

=  1^.,  y(p*(t))dt  . 


In  view  of  the  above  remarks. 


y„(P  ) 

T 


Y„(R)  >  P  (T-t)  -  P  (T) 
T 


(y(P  ( t)  )  -  si  dt  t  A 


* 

and,  since  P  (t)  +  P  , 

S 

lim  inf  (Y  (P*)  -  Y  (R)l  >  A  >  0 
T  T 

T  ♦  «*> 


(4h) 


For  A  >  0  this  result  and  Assertion  4  imply  (38),  Thus  assume  A  *■  0.  Since  R  ^  P* 
there  is  at  least  one  interval  Q  in  I  with  R  constant  at  S#  so  y(R)  -  ER  * 

S  -  ER  >  0  on  Since  P  =*  P  6  A  on  I  and  R  fails  to  satisfy  (A3)  on  c  i, 

P  /  R  on  il.  Thus  (38)  follows  from  (46)  and  Assertion  4,  This  completes  the  proof, 
d •  Discounting . 

The  usual  method  of  treating  an  infinite  horizon  problem  is  to  introduce  a  discount 
factor  e  into  the  yield.  The  total  yield 

f"  e"’^E(t)P(t)dt 

is  then  defined  for  all  bounded  paths  P  and  may  be  written  as 

Y(P)  =  /“  e"’^w(P(t))dt  +  Pj, 


where 


w(P)  *  y(P)  -  rP  .  (47) 

The  corresponding  optimal  harvesting  problem  for  unconstrained  effort  consists  in 
maximizing  Y  over  the  class  of  all  bounded  paths  which  arc  consistent  with  (A1)  and 
fA2),  If  w  is  a  function  on  [0,<»)  with  a  strict  global  maximum  at  >  0  and  if 

w  has  no  other  extrema  on  (0,»)  (48) 

then  the  solution  of  this  problem  is  analogous  to  the  solution  of  the  nondiscounted  problem 
of  section  b.  (cf.,  e.g.,  Clark  and  Munro  (1)  and  Clark  (2],  chapter  2.)  We  now  show,  by 
example,  that  if  the  rather  restrictive  assumption  (48)  is  eliminated,  the  optimal  path 
ne#=»d  not  be  a  most  rapid  approach  to 
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P  P  e 

0  0 


Figure  15.  The  function  w. 


Let  Pjj  >  0,  <  >  0,  and  M  >  0  be  fixed.  Further,  let  e  >  0  be  arbitrarily  small, 
let  “  "e  ®  function  on  [O,*)  with  the  following 

properties:  w  has  a  strict  global  maximum  at  P|^,  w  >  0  on  [Pq,P|^1, 

w(P„)  »  M,  w(Pq)  -  M  -  e  ,  (49) 


M  -  e  on  [P(j>Pge  1 


1  2  2 
€ 

e  on  (P  e  ,  P  e  1 

U  n 


we  will  show  that  for  e  sufficiently  small  the  path  P  =  is  better  than  the  path  P 

defined  in  ( 7 ) . 

To  see  this  let  Z  =  Zp  be  the  free-growth  curve  starting  from  P^  at  t  -  0,  so 


that ,  by  ( 47 ) , 


Z  =  w(Z)  ♦  icZ  . 


Let  t^  =  time  at  which  Z  reaches  P^^  (cf.  (7)).  Further,  let  tQ  *  ^0^®^ 

p2  _p2 

and  t^  =  t^(@)  denote  the  times  at  which  Z  reaches  PqS  and  P^^e  ,  respectively, 

so  that  t^  <  t^  <  t|^.  Assume  that  @  <  tc.  Then  (51)  yields 


PZ  i  Z  <  M  4-  ir-Z  on  r  n ,  t^  1 


-@2 

r-inre  M  and  <  are  Independent  of  e,  while  P^^e  becomes  infinite  as  S  *  0,  we 
nay  conclude  from  the  second  inequality  in  (52)  that 

t^+ooase+O.  ("'3) 

e2 

Similarly,  the  first  Inequality  in  (52)  yields  Z(tg)  >  P^e  ,  and,  since  Z(to)  “  ^0®  ' 

it  ‘’allows  that  t  <6;  thus 
0 


tjj  0  as  e  ♦  0  .  (54) 

In  view  of  (50),  w(Z{t))  is  kxjunded  above  by  M  -  e  for  0  <  t  <  t^,  by  e  for 

t  <  t  <  t  ,  and  by  M  for  t  <  t  <  t  >  thus,  using  (7)  and  (49),, 

0  1  1  M  ^ 


•JCt, 


y(P*)  -  Y(P)  <  /  ^  (2e-M)e“‘^dt  +  /"  ee''^dt 
^0 


(2e-M)e 


+  (M-3e)e 


By  (53)  and  (54)  the  right  side  of  the  above  relation  has  the  limit  -M/<  as  e  ♦  0, 

A  * 

Hence  we  have  only  to  choose  e  sufficiently  small  to  insure  y(P)  >  y(P  ). 


3.  The  age-dependent  theory  without  harvesting, 
a .  Basic  equations.  Persistent  age  distributions. 

We  work  within  the  framework  of  the  nonlinear  theory  developed  in  [5).  This  theory  is 
based  on  the  equations 


da  dt 

P(  t)  =  /”  p(a,t)da  , 


B{t)  =  p(0,t)  *  B(a,P(t) )  p(a,t)da  , 


(55) 


where 


c(a,t)  is  the  age  dlstrthutton  (the  number  of  individuals,  per  unit  age, 
of  age  a  >  0  at  time  t  >  0)i 
P(t)  is  the  total  population; 

B(t)  is  the  birth-rate? 
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y  is  the  death  function  (y(a,P)da  is  the  probability  of  dying  in 


the  age  Interval  (a,a+da)  when  the  population  is  P); 

0  is  the  birth  function  (0(a,P)  is  the  expected  number  of  births 
to  a  person  of  age  a,  per  unit  time,  when  the  population  is  P) . 

Here,  however,  we  assume  that  dependence  on  total  population  is  confined  to  the  death 
function.  More  specifically,  we  assume  that^ 


B(a,P)  B(a)  , 

u(a,P)  »  u  (a)  +  u  (P) 
n  e 


(56) 


In  (56)^,  u  (a>da  represents  the  probability  of  dying  of  natural  causes  during  (arS+da), 
^  n  ■ 

while  u  {P)da  is  the  probability  of  death  due  to  environmental  factors  (crowding,  etc.) 
e  - 

during  the  same  interval.  In  view  of  this  interpretation, 


T!(a)  =  exp{'  Vi^(a)da} 

is  the  probability  of  living  to  age  a  in  a  perfect  environment?  that  is,  an  environment 
with 

-  0  .  (57) 

e 

Assumption  (56)2  is  special,  as  the  effects  of  the  environment  are  independent  of  age?  it 

might  be  appropriate,  for  example,  to  a  population  in  the  presence  of  predators  which 

indiscrlminatly  eat  individuals  of  all  ages. 

We  assume  u  ,  y  ,  and  0  are  functions  on  [0,»)  with  B  H  and 

n  e 

/q  n(a) B(a)da  >  1  .  (5P) 


1 

These  assumptions  with  =  0  and  B(a)  a  sum  of  terms  of  the  form  were 

\itillzed  by  Gurtin  and  MacCamy  (5,8]  to  reduce  the  system  (55)  to  ordinary  <1i  f  f  erent  ia  I 
equations.  The  more  general  form  (56)  was  introduced  by  Coleman  and  Simmes  [4]. 
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The  left  side  of  (58)  represents  the  net  reproduction  r^tte^  in  a  perfect  environment.  A 
consequence  of  (58)  is  that  in  such  an  environment  the  population  ultimately  grows.  This 
assumption  therefore  yields  a  population  which  is  amenable  to  harvesting. 

Because  of  (58),  the  equation 


/”  Tr(a)6(a)e“*'®da  =  1 


has  exactly  one  real  solution  r,  and  r  >  0.  We  call  r  the  natural  growth  rate.  For  a 
perfect  environment  the  function 

p(a,t)  =  CTi(a)e  *^*P(t)  (59) 


with 


P(t)  =  P(0)e'^^ 


and  C  chosen  so  that 


C  /”  it(a)e  *  t  ($0) 

is  a  solution  of  (55)  ((56)  and  (57)).  Me  call  such  solutions  persistent  age 
distributions.  Their  importance  is  that  they  represent  the  asymptotic  behavior  of  general 

2 

solutions  for  large  time:  given  any  solution  p  there  is  a  constant  Cq  such  that,  for 
each  a , 

-ra  rt 

p(a,t)  ^  C^iT(a)e  e 

? 

ns  t  +  ». 


1 

Cf.,  Kpyfitz  [9],  p,  102. 

2 

1 1  er  [  1  . 

3 

f(^l  fj(t)  iq  t  ►  «  signifies  that  f(t)/q(t)  -*-1  as  t 
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Once  B  is  known  the  total  population 


P(t)  -  /”  C(a,t)da 


(65) 


corresponding  to  5  is  easily  calculated  using  the  relations 

P(t)  =  /J  «(a)8(t-a)da  +  'i'(t)  , 

f(t)  =  /”  ii(a-t,a)v>  (a-t)  da  . 


(66) 


It  is  important  to  emphasize  that  5,  8,  and  P  are  the  age  distribution,  birth-rate,  and 
total  population  that  would  prevail  in  a  perfect  environment. 

To  derive  (64)  and  (66)  from  (63)  we  integrate  (63)^  along  characteristics  (t  =  a  + 
constant).  The  result. 


5(a.t)  = 


{ir(a)8(t-a),  0  <  a  <  t 

III  a-t,a)'P  (  a-t)  t  6  a 


(67) 


when  combined  with  (63)2  and  (65),  yields  (64)  and  (66). 

Mext,  we  Integrate  (62)  with  respect  to  age  and  use  (55)2  and  (65)  to  conclude  that 


P(t)  =  P(t)exp{/p  yg(P(X))dX} 


(68) 


If  we  differentiate  this  relation  with  respect  to  t  and  divide  the  result  by  P(t),  we 

find  that  P  satisfies  the  ordinary  differential  equation 

P  =  (T  -  U  (P)1P  .  (69) 

e 


^(t) 

P(t) 


Here 

y(t) 

represents  the  instantaneous  growth  rate  in  a  perfect  environment.  Of  course,  P 
satisfies  the  initial  condition 


(70  ) 
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f*v(«)da  . 


(71) 


P(0) 


The  above  analysis  yields  the  following  procedure  for  solving  the  initial-value 
problem  (55),  (56),  and  (61); 

(a)  Solve  the  linear  integral  equation  (64)  for  B(t). 

(b)  Compute  P(t)  from  (66)  and  y(t)  from  (70). 

(c)  Solve  the  differential  equation  (69)  -  with  initial  condition  (71)  -  for  P(t). 

(d)  Compute  5(a,t)  from  (67)  and  p(a,t)  from  (62). 

4.  Harvesting  of  age-structured  populations, 
a.  Basic  equations.  The  optimal  harvesting  problem. 

We  consider  harvesting  with  effort  E(t)  independent  of  age.'  Thus  individuals  of 
age  a  are  harvested  at  a  rate  E(t)p(a,t),  per  unit  age  and  time,  and,  in  place  of 
(55)^,  we  have  the  balance  equation 

'^3'a~'*^~  ~^3t~^^  *  !u(a,P(t))  +  E(t))o(a,t)  =  0  .  (72) 

The  underlying  equations  are  therefore  (72)  and  (55)2  3  supplemented  by  the  initial 
condition  (61)*  (We  continue  to  assume  that  u  and  satisfy  (56)t) 

We  will  use  the  procedure  discussed  in  Section  3  to  reformulate  the  problem.  We 
therefore  define 


f:(a,t)  =  o(a,t)exp{[J  +  E(l>)dX)  .  (71) 

Then  ^  satisfies  the  linear  system  (63),  fv(t)  =  ^(0,t)  is  again  the  solution  of  the 
linear  integral  equation  (64),  and  r(t),  defined  by  (65),  is  again  given  by  (66).  It  is 


1 

Cf.  Remark  3. 
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important  to  note  that  B,  and  P  represent  the  age  distribution,  birth-rate,  and 


total  population  that  would  prevail  in  the  absence  of  harvesting  and  environmental  effects. 

As  before,  the  behavior  of  the  actual  population  P(t)  is  governed  by  an  ordinary 
differential  equation.  Indeed,  the  exact  same  steps  used  to  derive  (69)  now  lead  to 

P  =  y(t,P)  -  EP  (74) 

with 

y(t,P)  =  (y(t)  -  p  (P)]P  .  (75) 

e 

In  view  of  (74),  the  yield 

/g  /“  E(t)p(a,t)dadt  -  E(t)P(t)dt  (76) 


on  -i  time  interval  10,T)  is  given  by  the  functional 

=  Pq  “  ^  y(t,P(t))dt  ,  (77) 

Koughly  speaking,  our  optimal  harvesting  problem  consists  in  maximizing  (77)  subjet  to 
certain  constraints.  This  problem  can  be  attacked  using  the  following  procedure: 

(a)  Solve  the  Linear  integral  equation  (64)  for  B(t). 

(b)  Compute  r(t)  from  (66)  and  Y(t)  from  (70), 

(c)  Solve  the  optimal  control  problem  defined  by  the  differential  equation  (74) 
and  the  fvinctional  (76), 

nnfort’.xnatel Y,  because  of  the  dependence  of  y(t,P)  on  t,  the  optimization  problem 
of  (c)  ioe-^  not  fit  within  the  framework  discussed  in  Section  1,  There  is,  however,  an 
important  class  of  problems  for  which  this  dependence  disappears.  This  class  corresponds 
to  initially  persistent  age  distributions^  and  we  shall  study  it  in  the  next  section, 
h.  optimal  harvesting  for  initially  persistent  age  distributions. 

As  noto(i  at  the  end  of  Section  3a,  when  the  environment  is  perfect  and  harvesting 
,  f-h*'  3'ie  dl.9tr  ibut  ion  is  ultimately  persistent.  Therefore,  if  the  population  has 
been  Hvnlvina  ov(?r  a  Iona  period  of  time,  and  If  all  past  harvesting  has  been  with  effort 
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independent  of  age,  then  it  seems  reasonable  to  assume^  that  Initially  p  will  have  the 
age  structure  indicated  in  (59).  We  therefore  assume  that  the  population  Is  Initially 
persistent  in  the  sense  that 

■P  (a)  »  CPgV(a)e  , 

where  C  and  Pg  are  constants  with  C  chosen  so  that  (60)  holds,  while  r  Is  the 
natural  growth  rate.  The  Initial-value  problem  (63)  then  has  the  unique  solution 


5{a,t)  -  Cit(a)e“’^*Rt) 
Pit)  -  P^e*^*^  . 


(78) 


so  that,  by  (70), 


Y(t)  -  r  .  (79) 

Further,  the  counterpart  of  (68)  In  the  present  circumstances  Is  simply  (68)  with  E( \) 

added  to  y^(P(X));  thus  (78)^  and  (73)  Imply  that 

p(a,t)  »  Cit(a)e””p(t)  .  (80) 

In  view  of  (75)  and  (79),  y(t,P)  la  Independent  of  t: 

y(P)  -  (r  -  g  (P))p  .  (81) 

e 

Thus  the  differential  equation  (74)  reduces  to 

P  -  y(P)  -  EP  ,  (82) 


while  the  functional  (77)  takes  the  form 


“  *"0  ~  y(P(t))dt 


(83) 


Comparing  (82)  and  (83)  with  (3)  and  (4),  we  see  that  for  an  initially  persistent 
population  the  optimal  harvesting  problem  reduces  to  the  (age-independent)  problem 
discussed  in  Section  2.  Therefore,  if  y  satisfies  the  hypotheses  of  Section  2,  the 
theorems  of  Sections  2b  and  2c  hold. 


1 

Cf.  (73). 
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Remark  3.  Our  results  extend  trivially  to  situations  in  which  Individuals  of  age  a 


have  an  economic  value  g(a),  where  g  is  a  nonnegative,  L*  function  on  [0,“>).  In 
this  case  the  yield  (76)  is  replaced  by 


/q  /"  g(a)E(t>p(a,t)dadt 


/T  E(t)G(t)dt 


(84) 


where 


G(t)  -  /*  g(a) p(a,t)da  , 


and  (80)  implies  that 


G(t)  =  (C  /"  g(a)it(a)e  ’^*da)P(t)  =■  C^P(t) 

with  C,  >  0  constant.  Thus  (84)  is  a  constant  times  the  yield  (76),  and  the  corresponding 
optimization  problem  reduces  to  that  already  discussed.  When  the  age  distribution  is  not 
Initially  persistent  this  reduction  does  not  taVie  place;  the  resulting  problem  is  treated 
in  (11). 
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